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Abstract. In a previous paper [Phys. Rev. E 83, 051911] we have shown that the 
results of a quantum-mechanical calculation of electronic energy transfer (EET) over 
aggregates of coupled monomers can be described also by a model of interacting 
classical electric dipoles in a weak-coupling approximation, which we referred to as 
the realistic coupling approximation (RCA). The method was illustrated by EET 
on a simple linear chain of molecules and also by energy transfer on the Fenna- 
Matthews- Olson (FMO) complex relevant for photosynthesis. The study was limited to 
electronic degrees of freedom since this is the origin of coherent EET in the quantum 
case. Nevertheless, more realistic models of EET require the inclusion of the de- 
cohering effects of coupling to an environment, when the molecular aggregate becomes 
an open quantum system. Here we consider the quantum description of EET on 
a linear chain and on the FMO complex, incorporating environment coupling and 
construct the classical version of the same systems in the density matrix formalism. 
The close agreement of the exact quantum and exact classical results in the RCA is 
demonstrated and justified analytically. This lends further support to the conclusion 
that the coherence properties of EET in the FMO complex is evident at the classical 
level and should not be ascribed as solely due to quantum effects. 



PACS numbers: 82.20.Nk,82.20.Rp 
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1. Introduction 

In a previous communication [1] (to be denoted as paper 1) we studied theoretically 
the process of electronic energy transfer (EET) on molecular aggregates. The aim 
was to demonstrate that the coherent transfer arising from an entangled aggregate 
wavefunction (an exciton) in the quantum case, in the approximation that there is no 
exchange of electrons between monomers, is reproduced by a classical model of the 
aggregate as an assembly of electrical dipoles. The quantum/classical equivalence is 
valid in what we called the "realistic coupling approximation" (RCA). This is a weak- 
coupling approximation in that the strength of the dipole-dipole interaction that effects 
the transfer is considered small compared to typical electronic excitation energies, so that 
the monomers largely retain their character upon excitation of the aggregate. Practically 
this implies that the exciton bandwidth and the average spread of monomeric transition 
energies are both small compared to the mean electronic transition energy. Happily 
these criteria do pertain in many dye aggregates and also in the photosynthetic unit, so 
that the RCA is valid and one expects the classical model to give results in agreement 
with the quantum theory. 

As specific example, we considered first the transfer of energy along a chain of 
identical monomers where, in the approximation that only nearest-neighbors interact, 
an analytical solution for the transfer probability is possible. This solution predicts 
oscillatory monomer excitation probability in time and a constant velocity EET 
along the chain from an initially-excited monomer. The constant velocity and the 
oscillating nature of the transfer probability are signatures of fully- coherent propagation. 
Significantly it was shown, by numerical solution of the full classical equations for 
the same coupling, that classical interacting dipoles lead to exactly the same coherent 
transfer as in the quantum case. As a second model we considered EET on the Fenna- 
Matthews- Olson (FMO) photosynthetic complex, where the local transition energies of 
bacterio chlorophyll molecules on different sites are unequal. Here again, using realistic 
transition energies and coupling strengths, we were able to show the equivalence of 
classical and quantum dynamics in deciding the coherence of EET. The equivalence in 
this example is of particular significance since here the coherence of EET has been 
attributed [2] as arising solely from the entanglement properties of the aggregate 
electronic wavefunction. Hence, were this to be true, one would not expect such 
coherence to be present in the results of a purely classical treatment. 

In the above examples only electronic degrees of freedom were considered, since 
the emphasis was on coherent EET and the electronic excitation is the seat of such 
coherence. Nevertheless, particularly with reference to the FMO complex, in any real 
molecular aggregate the internal electronic degrees of freedom experience interaction 
with internal and external vibrational modes and electromagnetic interaction with the 
surrounding solvent. In absorption and emission of photons such interactions manifest 
themselves obviously in broadening and shifting of spectral bands. In EET there are 
more subtle manifestations in the de-phasing, de-cohering and sometimes dissipating 
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effects on propagation of absorbed light energy. Hence, to obtain a more realistic overall 
picture and to further test the classical model of EET we feel it essential to include 
the interactions with the surroundings. This enlargement of the theoretical model, to 
consider the molecular aggregate as an open quantum or classical system, is the subject 
of this paper. 

The development of the paper proceeds as follows. In section [2] we consider the 
quantum case of a molecular aggregate in which the coupling to the environment is 
represented by interaction of electronic degrees of freedom with external stochastic 
fluctuations. These lead to de-phasing and time-varying electronic transition energies. 
Here we adopt the usual density matrix formulation leading to a Lindblad-type master 
equation which is equivalent to that derived by Haken, Reineker and Strobl [3l H] 
arising from Markovian environment fluctuations. This will be denoted as the HRS 
equation. This model (and its variants) has been used extensively to describe exciton 
transport in molecular crystals, molecular aggregates and photosynthetic complexes (see 
e.g. Refs. [SIEIITIIHIEIITO]) 

In section [3] we apply the same physical assumptions to derive a classical density 
matrix equation of a similar, but not exactly equivalent, form to the quantum equation. 
In an appendix, the equivalence of the classical formulation in the RCA to the quantum 
version is proved. The proof is most transparent using the stochastic Schrodinger 
quantum equation (which is equivalent to the HRS equation). 

In section m the two cases already considered in paper 1 as 'bare' electronic systems 
are re-calculated including the effects of environment coupling. The main effect of this 
coupling is to damp out oscillations in occupation probabilities and to slow the rate 
of EET. Perhaps more importantly, by plotting density matrix elements we show that 
coherences between different sites are suppressed also. In the case of the FMO complex, 
our model can now be considered a realistic representation of the main features of EET 
and yet, significantly, the exact purely classical model including environmental effects 
again gives results hardly distinguishable from the fully quantum results, including the 
coherences between different sites. 

2. Quantum Mechanics 

The excitonic part of the aggregate is described by the Hamiltonian Hex = Hq + V 
where Hq is the sum of the Hamiltonians of non-interacting monomers and V is the 
total potential energy of the pairwise interactions between monomers. Since we consider 
the propagation of a single electronic excitation along the aggregate, we expand the 
Hamiltonian with respect to states | vr„ ) in which monomer n is electronically excited 
and all other monomers are in their ground state. In this basis one has 

Hex = '^(^n\'"'n){T^n \ + Km| T^n ) ( | (l) 

n n,m 

where e„ is the single- monomer transition energy and the full aggregate ground-state 
energy is set to zero. The matrix-element Vnm describes excitation transfer between site 
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n and m. 

In the following we are interested in the dynamics of the (reduced) density 
matrix p{t) of the electronic system when the electronic excitation interacts with an 
environment. We adopt a particular simple model where the dynamics of the density 
matrix p{t) is determined by a Lindblad Master equation of the form 

Pnm{t) = V.[p]nm + C[p]nm- (2) 

where 

'H[p]nm= - ^[Hc^,p{t)]nm (3) 

= — ^i^n — ^m)Pnm — ^^^(Ki^Pfm — yimPni) (4) 

e 

and the last line follows from ([1]). The interaction with the environment is contained in 
^[p]nm which for simplicity we take to be 

'Clpjnm = - (^(7n + 7m) " V^nlmSnm^ Pnm (5) 

i.e. we consider pure dephasing with dephasing rates 7„. 

The extension to the general HRS master equation is straightforward. Also 
generalizations of the HRS model, as in Ref. [HI [121 CHI E], can be treated similarly. 
In the Conclusion we will discuss this point in more detail. 

In the following we will derive a classical equation which in RCA is equivalent to 
Eq. (|2]). To this end we first note that the master equation ([2]) is equivalent (following the 
treatment of HRS [3111]) to a stochastic Schrodinger equation dt\ il^{t) ) = —iH(t) \ ip(t) ) 
with local Markovian fluctuations of the site energies. These fluctuations can be merged 
together with the transition energy of the monomers to obtain a stochastic Hamiltonian 

H{t) = ^ e„(t) 1 7r„ ) ( 7r„ I + ^ Km| 7r„, ) ( 7r„ | (6) 

n n,m 

where the fluctuations in the transition energies e„(t) have the properties of real Gaussian 
Markov processes fulfilling 

iUt)))=en (7) 

(( en{t)em{t') )) = h\,6nm6it - t') (8) 

where ((■■■)) denotes the averaging over many realizations of the stochastic processes. 
The master equation ([2]) is then obtained by taking the time derivative of p{t) = 
{{\4'(t)){ip(t) \)). We note that we have restricted ourselves to Gaussian Markov 
processes, since then one can derive a simple master equation. The arguments presented 
below on the validity of the RCA approximation are applicable to more general stochastic 
processes which do not have to be Gaussian (as long as the second moment exists) and 
can also be correlated in time (non-Markovian). 

The stochastic 'unravelling' (I6|) will be used in the next section to make the 
connection to the classical case. In particular we will take the frequency of the classical 
oscillators to obey the same statistical properties as the quantum transition energies. 
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The resulting classical stochastic equation will then be used to construct a 'classical 
master equation'. 

3. Classical mechanics 

As in our previous publications [H [15] in the classical case we consider that the coupled 
quantum transition dipoles are modelled by classical oscillators in the same geometry 
as the transition dipoles of the quantum mechanical model. The frequencies of the 
classical oscillators are associated with the quantum energies via Un = ^n/h. To arrive 
at a density matrix description corresponding to ([2]) we make use of the stochastic 
representation introduced in the previous section. Thus we assume that the frequency 
of the classical oscillators is a stochastic quantity with 

UJnit) = en{t)/h (9) 

and the same statistical properties as given by (I7j) and ([H]). 

The Hamilton equations of motion for linearly-interacting oscillators of time- 
dependent fiequencj Unit) as defined above in ([9]), are, 

Xn = UJn{t)Pn (10) 

Pn= - (^n{t)Xn - ^ KnmXni (H) 
m 

where x„ and Pn are the dimensionless position and momentum of the nth oscillator 
respectively. The Knm are coupling coefficients which are related to the quantum 
mechanical couplings by [H [15] 

= 2^. (12) 

To connect to the quantum equations, we introduce the dimensionless complex 
amplitude 

Zn{t) = Xn{t) + ipn{t) (13) 

and obtain the coupled equations 

Zn = -iun{t)zn - i ^ ^^2Re(^m) (14) 

m 

These equations can be viewed as a set of coupled "Kubo-oscillators" [161 HI]- Note 
that here 2Re(2;m) appears as coupling in the equation for the amphtudes. As shown 
in [1] in an equation which would be fully equivalent to quantum mechanics this term 
would be replaced by Zm- 

3.1. Classical density operator 

To make contact with the pure dephasing master equation ^ we will take a closer look 
at the products 

C^nm{t) = Zn(t)z^(t) (15) 
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which, as we will show below (in RCA) resemble the quantum mechanical density 
operator matrix elements. In the following derivation some care has to be taken due to 
the stochastic nature of Un(t). We interpret (fT4l) to be a stochastic Schrodinger equation 
in the Stratonowich form (see e.g. Refs. [HI [T71 [IS]). In the following derivation we use 
the Ito calculus (see e.g. Ref. [18]) and write (HM in its Ito form 

dZn= iUnZn - i ^ Knm^eZm. - y^n) + dWn (16) 

m 

Here dz„ is the increment of during the time interval dt and dWn are Wiener 
increments fulfilling (( dVF„dVFm )) = ^nmdt and (( dWn )) = 0. Note that the first term 
of (ITB]) the oscillator frequency aj„ is the average frequency and does not depend on 
time. The effect of the stochastic fluctuations are contained in the factors ^Zndt and 

y/%ZndWn. 

We now derive the equation of motion for a. To this end we consider the differential 

of (USD 

danm = dznz*^ + z^^z*^, + dz^dz*^ (17) 

= iUnZn - i ^ knm'^eZm' - y^n.jdt Z*^^ + y/%:,ZnZ*^dWn 



+ Zn{iUJmZ*m + ^ Kmm''^eZm> " '^zl^dt + ^/^ZnZl^dW J^l%) 



+ ^lnlmZnZ*^dWndWm- 

Here we have taken terms up to the first order in dt into account. Since dW scales like 
\/di in the first line we have included the dzndz"^. 

We are interested in quantities that are averaged over the noise, where, as before, 
the averaging is denoted by ((■■■))• Defining the classical density matrix 

O-nm = {{^nm)) (19) 

we find, using ([7]) 

danm = - ii(^n " w„) " ( y + y))o-nmdt 

Re.™' z:, ))dt + ^ J2 ^(( ^nRez^, ))dt (20) 

m' m' 
^~ y/'^n'ym^nm^nmdt 

The term in the last line results from the averaging of the expression containing 
dWndWrn. 

If one compares (l20l) with (jlj) one sees that the the term —i{uJn—uJm)o'nm corresponds 
to —{i/h){en — €m)pnm- The tcrms containing 7's can be combined to give C[a] where 
C is the same functional as in the quantum case. The remaining terms contain the real 
part of the complex amplitude and therefore cannot be written as o"„m- We will now 
first bring (l20l) into a form which is closer to the quantum equation and then show that 
in RCA they become identical. Using 2 ReZm' = ^m' + z^, we can re- write (j20l) to obtain 

d-nm = 'H[0-]nm + C[a]nm + i ^ (( ZiZn )) - (( ))) (21) 
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This equation has to be compared to the quantum mechanical master equation ([2]) . Note 
that fl2T]) as it stands is not a closed system of equations for a, due to the appearance 



of (( Z£Zn)) and {{z£Zm))*- In [Appendix A| it is shown how we solve this equation. 

We still have the freedom to normalize the classical density operator. This will be 
done by the factor 

Ar=J2(^nn (22) 

n 

SO that we can identify 

p{t)^a{t)/Af (23) 

With this we have related the classical master equation to the quantum master equation. 
In the following subsection we will briefly discuss the initial state. 

3.1.1. The initial state 

Pure states. Consider first a quantum mechanical initial state of the form 

p'"^^ = \^'^'){^'^'\ (24) 

Writing the initial wave-function as | -j/;'"^ ) = ^„ c™^| Tin ) with |cJ^'P = 1 we get for 
the matrix elements of the initial density operator 

Pnm ~ i^m ) (25) 

The corresponding classical initial state is constructed by choosing 

C = «C (26) 

where a is an overall constant which will drop out in the end, when calculating 
populations, coherences, etc. Thus we have for the elements of the initial classical 
density matrix 

= «'C(C)* (27) 

Mixed states. To treat mixed states we first note that an arbitrary density matrix can 
be written as a weighted sum of pure states 

r = Y,wp\^Pp){^p\ (28) 

with p| "^A/j ) = wplipiB)- This suggests to construct the corresponding initial classical 
state as 

^'" = E^/5^/5 (29) 
where (cr/j)^^ = a'^{pp)nm- 
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3.2. The realistic coupling approximation 

To investigate under which conditions the RCA will be valid we use similar arguments 
to those in our previous work [H [15]. To this end we will consider not the density 
matrix equations but the equivalent stochastic equations. Expanding the quantum wave 
function as 

l*W) = E^«WI^") (30) 

n 

in the quantum case one has from ([6]). 

c„ = -i(e„ + hn{t))cn - i ^ KmC^ (31) 

m 

which has to be compared to the classical equation (0311 which can be written as 

Zn = -ijuJn + Wn{t))Zn - i ^^^Re(Zm) (32) 

m 

where the stochastic processes can be chosen to be Wn{t) = hn{t)/h and as before 
In Ref. [1], where h{t) = 0, it was noted that under the conditions that 

\Vnm\/fi <W„ (33) 
|W„ -UJm\ < (34) 

the classical equations accurately describe the results obtained from the Schrodinger 
equation. Then also the corresponding classical master equation should give results 
which agree with those obtained from the quantum master equation. 

In the present case, where h{t) ^ 0, we expect that the difference between 
the quantum and the classical evolution is small if the same conditions for the time 
dependent frequencies are fuUfilled, namely that 

I Vnm \/h<^UJn + Wn{t) (35) 
Un \Un + Wn{t) - {Ujm + Wm{t))\ (36) 

for most times t. Since Wn{t) can take negative values we see from Eq. (l35l) that the 
fluctuations have to be small compared to to the average frequency cu^. The second 
equation (l36l) states that the frequency difference between different sites has to be small 
compared to the mean frequency. 

For the situation considered in the present work we use 7„ as a rough measure of 
the magnitude of the frequency variations. Then both equations ( l35l) and ( l36l) lead to 
the estimate 

In -C Un (37) 

If this condition and ( l33i) and ( l34l) hold, then we expect the RCA to be a good 
approximation. Note, that these assumptions have to be fuUfilled also for the quantum 
mechanical model employed to be meaningful. That is the energy changes experienced 
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by a given molecule due to coupling to other molecules and the environment must be 
small compared to the unperturbed molecule transition energy. 

In section m we will investigate the range of validity of this approximation in more 
detail. In particular we will show that for typical parameters used in the photosynthetic 
FMO complex the classical equations give a good description to the quantum mechanical 
exciton dynamics. 



4. Comparison of quantum and classical results 

Here we compare the results obtained using the classical "master equation" ( 12T|) with 
those obtained from the quantum one ([2]). The classical equation (12T1) is solved using 



the method described in [Appendix A 



First we will consider the case of a linear chain and afterwards discuss the FMO 
complex. For both cases we investigate the populations as well as the coherences between 
different sites. 



4.1. The linear chain 

As a first example we consider the standard case of a linear chain where all transition 
energies e„ and dephasing rates 7„ are identical. For simplicity we take only nearest 
neighbor interaction into account. This interaction, denoted by V is taken as the unit 
of energy. 

In the results shown in figure [H we have used as initial condition a state where the 
excitation is initially localized on a single monomer which we denote by '0' In each 
column a fixed 7 is chosen and the transition energy e is varied. For reference purposes, 
in the left column, the quantum solution for 7 = is also shown. Here the coherence 
is maximal and the probability of EET is given by the square of a Bessel function [20] 
which is oscillatory in time, reaching zero at the zeroes of the Bessel function. As 7 is 
increased one sees that the oscillations are damped for j = V (left column of Fig. [1]). 
When 7 = 20V (right column) the oscillations are washed out completely and the 
populations change monotonically in the quantum case. The exact quantum results are 
shown as solid black lines. 

To investigate the validity of the RCA we show classical solutions for various values 
of the transition energy e = hw. In particular we have chosen e = 40, 30, 10, 6, and 1 
in units of V. For the values e = 40, 30, 10 Eq. (l33l) is fulfilled while it definitely does 
not hold for e = 1. Note that for our choice of identical transition energies Eq. (IMIl 
is trivially fulfilled. For the case 7 = 1 we have 6/7 = e/V and the inequality ( 1371) is 
fulfilled whenever fl33|l is fulfilled. However for 7 = 20 the inequality fl37j) is not fulfilled 
for e = 10, 6, 1. For this case one expects to see the influence of the fluctuating transition 
energies on the RCA. 

I For the calculation wo used a chain of 29 sites and started at site 14. 
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Figure 1. Time-dependent populations of a linear chain when the excitation is 
localized initially on one monomer "0". The transport is symmetric w.r.t. this 
monomer. Left: 7 = y (for comparison also the analytic 7 = result is shown as thin, 
gray curve). Right: •y = 20 V. Bold black: exact quantum calculation. The colored 
curves are results from the classical calculations for different e. Orange: e = 40 V. 
Blue: e=10V. Red: e^6V. Green: e^lV. Time is in units of V/h. 



These expectations are met by the numerical results shown in Fig. [H We see that for 
the case 7 = 1 for e > 10 the quantum and classical results are nearly indistinguishable. 
Even for a ratio 5/7 = 6 (red curve) there is still quite good agreement. For even smaller 
e/7 the deviations become more pronounced as exemplified by e/7 = 1 (green curve). 
For the large value 7 = 20 the classical results are indistinguishable from the quantum 
ones for the cases e > 30 (where both inequalities ( 137|) and (l33l) hold). For all other 
values (where e < 7) there are clear deviations from the quantum result. 

Not only the populations but also the inter-site coherences obtained from the 
classical 'density matrix' are in good agreement with the quantum mechanical 
coherences. This is demonstrated in figure [2]for the case 7 = V where in the left column 
for a linear chain the time- dependent absolute values of the coherences between site '0' 
and site '1' (upper row), site '2' (middle row) and site '3' (lower row), obtained from 
full quantum calculations, are shown. Oscillations in the coherences are evident. Since 
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Figure 2. Left column: Quantum coherences |po,i|i |Po,2| and Ipo.sl of a linear chain 
when initially the excitation is localized on one monomer "0" . The dephasing rate is 
7 = y. Middle column: Differences between classical and quantum results for the case 
e/7 = 40. Right column: Differences between classical and quantum results for the 
case e/7 = 6. Note the different scalings in columns (b) and (c). Time is in units of 

v/h. 



deviations between classical and quantum results are not always easy to distinguish, in 
the middle column the differences between the absolute values of the quantum result 
and the classical calculation for the case e/7 = 40 are shown. One sees fast fluctuations 
in the differences but they are always at least two orders of magnitude smaller than the 
magnitude of the exact coherences. Upon decreasing the ratio e/7 the fluctuations in 
the difference become larger and slower but never exceed more than a few percent. An 
example is shown in the right column for e/7 = 6. 

4.2. The photo synthetic FMO complex 

The simple HRS model can be used to gain insight into the dynamics of excitation 
energy transfer in photysynthetic systems, although for a more realistic description a 
more detailed treatment of the environment would be necessary [211 1221 |23l [23] . 

Here we show, for the case of pure dephasing (with 7„, = 7), that for values of 
the transition energies and interactions between the BChl molecules that are typical for 
such systems the classical equation gives results in very good agreement with those of 
the quantum one. 

In the following we use the energies and couplings for c. tepidum as given in Ref. [25] 
(the site energies are taken from the trimeric structure of table 4 and the intersite 
couplings are taken from the fourth column of table 1 of that paper). We note that for 
our comparison the exact values are not important, therefore we did not consider more 
recent values [261 [23] . One would get similar results as presented below for slightly 
different parameters and also when treating the full trimeric FMO system with 24 BChl 
molecules. 
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Figure 3. a) Populations of the BChls as function of time when excitation is initially 
localised on BChl 1 obtained from the full quantum mechanical calculations, b) 
Differences between the quantum result and the classical calculation for the actual 
transition energies, c) The transition energies are reduced by 12000 cm~^ . Note the 
different scales of the y axis and the scaling factors which are indicated in the figures. 

In figure |3]the time dependence of the excitation probabilities of the Bchl molecules 
is shown. The left column shows the results from the full quantum calculation. The 
middle and the right columns show the differences between the quantum result and 
the solution of the classical equation. In the middle column the transition energies are 
taken as given in Ref. [25]. They are in the order of 12000 cm~^ which is much larger 
than the coupling between the BChls (which is in the order of 100 cm~^) and also much 
larger than the energy differences between the transition energies (which are also a few 
hundred wavenumbers). As expected the the deviations from the exact quantum result 
are quite small (on the order of 0.1%). Even if one reduces the transition energies by 
12000 cm~^ (which results in completely unrealistic transition energies of the order of a 
few hundred wavenumbers) the deviations are still only on the order of 10%. 

The agreement between the classical coherences and the quantum coherences is of 
a similar order. This is exemplified by the results shown in figure HI The absolute 
differences are of the same order of magnitude as the differences in the populations. 
However, since the coherences are smaller, the relative error is somewhat larger. 
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Figure 4. As figure [3] but now for tfie colierences of BClil 3 witli all the other BChls. 
Note that in the third row also the diagonal term i.e. the population of BChl 3 is 
shown. Note the scaling factors in column (b) and (c). 



5. Conclusions 



In the present paper we have extended our previous investigation on the correspondence 
between quantum mechanical and classical EET to include coupling to an environment. 
In particular we have demonstrated that it is possible to derive a master equation for the 
classical amplitudes which in the RCA reproduces the corresponding quantum master 
equation. This has been shown explicitly for the case of an environment that leads to 
pure dephasing. As shown by HRS the corresponding quantum master equation can be 
obtained from an average over a stochastic Schrodinger equation with real Markovian 
noise that alters the transition energies of the monomers. We used this stochastic 
representation to relate the quantum Schrodinger equation to a classical equation for 
coupled harmonic oscillators with frequencies that have the same stochastic properties 
as the quantum transition energies. We then found that the classical results reproduce 
the quantum results when the fluctuations in the transition energies are small compared 
to the transition energy. This has been demonstrated explicitly by considering both a 
linear chain and the FMO complex. 

Although we have made this demonstration for real Gaussian Markovian noise it 
is clear that the same argument will also hold for more general stochastic processes for 
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the transition energies, which may be non-Gaussian or non-Markovian (of course also 
deterministic functions and fluctuations of the couplings between different sites can be 
treated in this manner). 

In this respect it is worth mentioning that such fluctuating site transition 
energies arise in molecular dynamics/quantum chemistry simulations performed on 
photosynthetic complexes [271 [221 [2H]- Thus the transport and coherence properties 
obtained from such studies could be reproduced by a purely classical model. 

The results presented demonstrate that classical oscillators can be used to simulate 
the transport and coherence properties of molecular aggregates and in particular those of 
photosynthetic systems. In a recent article [I5] we have discussed how classical electrical 
LC circuit oscillators can be used to mimic coupled quantum two level systems. From 
the foregoing it is clear that by modulating the frequencies of the classical oscillators 
in an appropriate way the classical LC oscillators can be used to simulate the quantum 
aggregate. Such classical simulations could be used to compare with recently proposed 
simulations using superconducting qubits [29] . 
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Appendix A. Solution of the classical equation 

In this section we describe how we solve the classical equation (j2T]) . To this end we 
introduce the auxiliary matrices 

Snm = {{PnPm)) (A. 2) 

Tnm = {{S:nPm)) (^-3) 

With this the classical density operator can be written as 

From the evolution equations of x„ and pn we can then derive a set of coupled 
equations: 



Rnm ^nTmn ~l~ ^mTnm ~\~ C,\R)[^nm (-^•^) 

'Vnerp '^m 
-^J-lm H T' 

a n 



Snm — — ij^uTnm + ^mTmn) + C,[S]nm — — T~Tim H T — ^te) (A. 6) 



Tnm ^nSnm ^mRnm ^\T^nm. ^ ^ Rni (■^■''') 
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Appendix A.l. The initial state 

In the following we set for convenience the normalisation factor a, introduced in ( l26l) . 
to a = 1. For each /3 in ( l29l) we make the identification 2;^ = c^. Then we construct the 
components i?^, and from the equations (]A.l|) - lA73l This gives 



Y^'^P'^nm (A.8) 

/3 

— _|_ Q™' _|_ ^YT'"' _ T™M I'A Ot\ 

^nm ~r "-^nm ~r ''\'''mn -'"nrn) \-^-^) 

with = E/3^/3^L. 'S'X = E/j^/J-^L and = E/3^/3^nm- Thus the initial 
vector (i?™', 5'"^ T™'), needed for the propagation with (12T]) . can be written as a sum 
of the vectors S*^, T^). The linearity of (lA.5l) -( |A7n) then guarantees that a{t) can 
be obtained at later times. 

Appendix A. 2. Comparison with the qm- equation 

Similarly as for the classical equation one can also re-write the quantum equation, giving 
further insight into the RCA approximation. In the quantum case we define 

R^Z = ReCnReCm (A. 10) 

= Imcjmc^ (A.ll) 

TZ =ReCnlmCm (A. 12) 
With this the density operator can be written as 

Pnm = Rnm ~^ ^nm + ^(-^nm + '^mn) (A. 13) 

We can then derive the set of coupled equations: 

TZ = ^nSZ - ^mRZ. + >c[r^-]_ + + Yfsz) (A.ie) 

One sees that compared to the classical case the terms containing the interaction between 
the monomers are more symmetrically distributed among the individual equations. 
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